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I INTRODUCTION 


Tliis  technical  report  is  intended  to  describe  some  of  the  work  done  at  the  Naval  Surface  Weapons 
Center,  Dahlgren  Laboratory  in  FY-76  on  the  application  of  Geometric  Theory  of  Diffraction  (GTD) 
analysis  and  related  methods  to  shipboard  electromagnetic  problems.1  In  particular,  the  work  addressed 
the  search  for  methods  which  are  accurate  enough  to  give  useful  results,  yet  simple  enough  to  permit  solu- 
tion of  the  very  complex  analysis  problems  arising  in  the  shipboard  electromagnetic  environment.  As  a 
“testbed”  problem,  the  scattering  of  electromagnetic  radiation  from  a metallic  circular  cylinder  was  chosen. 
Tliis  problem  could,  for  example,  represent  the  scattering  of  the  radiation  from  a radar  antenna  by  a ship's 
mast  or  a whip  antenna.  The  problem  to  be  solved  could  for  instance  be  as  follows:  We  install  a radar 
antenna  in  location  A on  board  the  ship  and  a cylindrical  mast  in  location  B.  We  intend  to  transport 
hazardous  explosives  along  path  c.  Will  the  path  c be  in  an  hazardous  electromagnetic  field?  If  yes.  how 
can  the  geometry  of  A versus  B versus  c be  modified  so  as  to  eliminate  the  hazard?  In  order  to  obtain  a 
good  appreciation  of  such  problems,  it  is  important  to  remember  that  there  is  a hierarchy  of  contributions 
to  the  total  field  in  any  location,  which  consists  of  a superposition  oft  I ) direct  radiation,  (2)  radiation 
reflected  from  the  cylinder  and  (3)  radiation  diffracted  from  the  cylinder.  In  the  following  sections  an 
attempt  will  be  made  to  establish  the  order  of  magnitude  of  these  contributions. 

The  estimates  described  in  this  report  and  in  the  proposed  follow-up  to  this  work  arc  designed  to  be 
worst-case,  i.e..  high  estimates.  The  computations  are  designed  to  serve  as  a first,  simple,  rapid,  and  cost- 
effective  cull  for  a large  number  of  field  points.  In  certain  regions  this  estimate  will  result  in  high  values  for 
the  total  field.  If  these  values  are  unacceptable  from  the  HERO.  RADHAZ.  EMC,  or  operational  point  of 
view,  a second  cull  will  be  performed  for  the  relatively  small  number  of  field  points  determined  by  the 
first  cull.  The  second  cull  will  employ  the  accurate  methods  developed  or  under  development  by  the 
academic  community.  This  two-step  process  will  optimize  the  efficiency  of  the  computational  procedures. 

The  effect  of  the  blockage  by  the  cylinder  may  be  divided  according  to  the  field  regions  considered. 
These  may  be  classified  as  shadow,  half-shadow,  and  unobstructed  regions.  In  the  work  reported  here  we 
concentrated  on  the  latter  two  regions,  which  in  most  practical  problems  will  make  up  most  of  the  total 
field  space  considered.  The  ratio  between  the  extent  of  the  total  shadow  region  on  the  one  hand  and  the 
half-shadow  and  unobstructed  regions  on  the  other  is  typically  between  the  orders  of  magnitude  0.01  and 
0.1  Although  the  problems  arising  in  the  shadow  region  are  sometimes  very  important,  we  did  not 
address  them.  In  the  two-step  process  described  above  this  region  of  the  electromagnetic  field  would  be 
handled  in  the  second  cull. 

The  work  was  accompanied  by  a literature  search,  which  also  permitted  us  to  place  our  work  in  per- 
spective with  respect  to  a representative  cross  section  of  the  published  literature.  Our  testbed  problem  is 
simple,  compared  to  the  much  more  complex  geometries  encountered  in  a shipboard  environment.  Yet  the 
literature  contains  a very  large  number  of  papers  which  in  some  way  treat  a subproblem  of  the  above  simple 
problem.  The  incident  field  may  be  a plane,  cylindrical,  or  spherical  wave  corresponding  respectively  to  a 
point  source  at  "infinity.”  a linear  source,  and  a point  source  “nearby,”  where  the  words  “infinity”  and 
“nearby”  have  a meaning  only  if  related  to  another  length,  most  frequently  the  wavelength  and/or  the 


diameter  id  the  e\  limler  Sometimes  complex  somces,  such  .is  arrays  ol  line  sources  or  special  sources  such 
.is  dipoles  are  treated  The  incident  field  may  be  further  specialized  by  its  time  dependence.  CW  pulse,  or 
step  function  The  problem  may  he  treated  as  two-dimensional  or  three-dimensional,  as  scalar  or  vectorial. 

The  secondats  field  heated  may  he  the  reflected,  diffracted,  or  refracted  field,  neat  field,  or  lai  held,  the 
interest  mas  enter  on  blockage  ol  the  mamlohe  or  generation  by  scatter  ol  new  sidelobes.  backscattcr  or 
Instatic  scattering.  01  special  regions  of  the  field  such  as  caustics  or  transition  regions.  A lurther  source  ol 
differentiation  lies  In  the  mathematical  methods  used  for  the  solution:  integral  ecpiations  method,  modal 
solutions,  conformal  mapping.  Ol  particular  interest  aie  the  possible  approximations.  There,  those  approxi- 
mations w hich  permit  simple  physical  analogies,  such  as  geometrical  optics,  physical  optics,  geometric  theory 
ol  diffraction,  are  particularly  appealing  because  they  are  helped  by  intuition  based  on  simple,  well-known 
physical  laws  The  physical  properties  ol  the  scattering  body  (cylinder)  range  from  minute  to  zero  con- 
ductivity li  e.,  metal  to  dielectric)  and  include  complex  conductivity  (e.g..  lossy  dielectric).  Depending  on 
the  relative  direction  of  the  polarization  of  the  incident  field  and  the  cylinder  axis,  the  boundary  conditions 
may  be  "soil”  or  "hard."  or  an  “impedance”  boundary  condition  may  apply.  Several  materials  may  be  in- 
volved in  the  scattering  process,  such  as  in  scatter  from  dielectric  covered  metallic  cylinders,  or  plasma- 
coated  cylinders.  Some  published  papers  center  not  on  the  scattered  field  hut  on  an  intermediate  quantity 
such  as  the  current  induced  on  the  surface  ol  the  cylinder.  Others  concentrate  on  mathematical  properties 
of  functions  occurring  in  the  treatment  of  scattering  problems,  e.g..  cylinder  functions.  Airy  integrals,  l ock 
functions,  finally,  sometimes  the  problem  treated  in  a paper  may  be  more  general  in  one  respect  and  more 
restrictive  in  another  than  our  testbed-cylinder  scatter  problem.  For  instance,  the  scattering  ol  a plane  wave 
front  a cylinder  of  arbitrary  cross  section  and  axis  perpendicular  to  the  incident  plane  wave  would  fill  into 
this  category. 

The  above  description  of  the  scope  of  the  literature  relating,  more  or  less  closely,  to  the  testbed  problem 
is  by  no  means  exhaustive  but  may  serve  as  an  illustration  of  the  possible  number  of  parameters  which  may 
combine  in  many  different  ways  in  the  literature. 

Considering  the  large  number  of  parameters,  the  "sample  space"  ol  possible  different  papers  which 
would  applv  to  the  different  possible  problems  occurring  in  actual  situations  is  enormous.  Ol  course  only  a 
small  fraction  of  this  sample  space  has  been  reported  on  in  the  literature. 

I would  like  to  end  this  introduction  by  mentioning  an  unbalance  I observed  in  scanning  the  literature. 
There  is  a large  number  of  purely  theoretical  papers  and  a far  smaller  number  of  experimental  papers.  This 
unbalance  is  at  least  partly  due  to  the  very  considerable  advances  in  theory  made  possible  in  recent  years  by 
the  enormous  progress  in  computer  science  and  available  computer  facilities.  However,  the  need  for  the 
interplay  between  theory  and  measurement  is  now  as  necessary  as  always. 


II  ESTIMATE  OF  REFLECTED  VERSUS  DIFFRACTED  POWER 


The  general  t.TD  solution  for  scatter  from  a circular  cylinder  was  published  many  years  ago  in  a paper 
by  I B Keller  el  al  2 It  involved  fairly  time-consuming  mathematical  operations  including  infinite  series. 

Our  purpose  was  to  obtain  an  extremely  simplified  formulation  and  al  the  same  time  an  estimate  of 
the  price  in  accuracy  ol  this  simplification  We  regarded  this  as  an  area  of  things  which  should  be  done  but 
which  do  not  show  up  in  detail  in  the  literature  Tiie  CTD  solution  to  the  cylinder  diffraction  problem  is 
based  on  an  asymptotic  series  tor  the  field  Cot  the  Luneberg-Kline  torm 


U(x,  e.  :)  = <>'* v 


whore  I m i\  thew-tb  order  term  ol  the  series,  and  k = 2 7r/X.  In  the  high-frequency  limit  the  series  is 
approximated  hv  the  zero-order  term  lor  m = 0.  Ii  depends  no  longer  on  ka  (where  u is  the  radius  of  the 
cylinder)  and  is  recognized  as  the  solution  ol  the  I’M  problem  offered  by  geometrical  optics.  As  ka  decreases, 
the  relative  values  of  the  higher-order  terms,  representing  diffraction  rather  than  reflection,  increase.  An 
examination  ol  the  expressions  lor  the  \'m  coefficients  shows  that  the  zero-order  lent!  is  dominant  in  our 
problem  for  most  held  points,  and  for  many  shipboard  environment  problems  only  the  first  term  in  the 
series  needs  to  be  considered.  This  is  true  for  both  normal  and  oblique  incidence,  l or  a plane  wave 
impinging  normally  on  a circular  cylinder  we  estimated  the  ratio  /?  between  total  diffracted  and  reflected 
energy.  We  use  l ock's3  expression  lor  the  width  d of  the  shadow  boundary  strip 

d = 

n 


where  a is  the  radius  of  curvature  of  the  body  in  the  plane  of  incidence:  i.e.,  in  our  case  a = radius  of  the 
ey  Under.  Assuming  that  the  diffracted  energy  is  due  to  the  energy  falling  onto  the  “region  of  the  shadow 
boundary"  as  defined  by  l ock’s  equation  and  assuming  further  that  the  optical  contribution  to  the 
scattered  energy  is  due  to  the  energy  tailing  onto  tne  part  of  the  cylinder  surface  lacing  the  arriving  plane 
wave,  we  may  arrive  at  an  estimate  of  the  relative  importance  of  the  total  diffracted  I.M  energy  to  the  total 
reflected  energy: 


defining 


reflected  energy 
diffracted  energy 


We  obtain  (see  Appendix  Al,  Equation  A 1 -3 ) 


R i = 


(i) 


where  R , is  an  approximation  for  the  desired  ratio.  Clearly  if/?,  is  large  we  shall  be  able  to  simplify  many 
problems  considerably  by  neglecting  the  higher-order  diffraction  terms  in  comparison  with  the  reflection 
term. 


Equation  (I ) is  evaluated  in  Appendix  Al  in  Table  Al  -1  for  various  values  of  ka  ( a is  the  radius  of  the 
cylinder)  between  ka  = 50  and  ka  = 2.  As  expected,  when  the  diameter  of  the  cylinder  decreases  or  the 
wavelength  increases,  the  relative  contribution  of  diffraction  in  the  total  scattering  process  rises  as  compared 
with  reflection.  However,  it  is  obvious  that,  since  the  thinner  cy’mder  interacts  less  than  a large-diameter 
cylinder  with  the  incoming  radiation,  the  total  diffracted  power  dt  .reases  with  cylinder  diameter.  This 
situation  is  shown  quantitatively  in  Appendix  Al , Table  A 1 -2.  Note  that  as  long  as  the  GTD  is  valid  (ka  > 2) 
the  total  reflected  power  dominates  the  total  diffracted  power. 


Although  the  analysis  in  Appendix  Al  assumes  a plane  wave  with  incidence  perpendicular  to  the  axis 
ol  a circular  cylinder,  it  can  be  extended  to  oblique  incidence  and  arbitrary  smooth,  convex  cylinder  cross 
sections  following  the  method  of  Appendix  A I and  substituting  for  the  radius  a the  radius  of  curvature  of 
the  cylinder  surface  in  the  plane  of  incidence  of  the  respective  points  in  the  shadow  boundary:  the  width  w 
of  the  penumbra  region  will  then  be 


tv  = 


Since  most  wavefronts  may  be  represented  by  a “plane  wave  spectrum.”  the  solutions  obtained  will  he 
applicable  to  a fairly  general  range  of  problems  involving  smooth,  convex  cylinders  in  near  or  far  fields. 


where  a radius  of  the  cylinder 

p.  cv  . . polar  coordinates  of  a point  ( on  the  caustic 
0.p  . ..  see  figure  III- 1 

Equation  (III  I I is  graphed  in  Figure  1111  above 

Note  that  the  virtual  caustic  is  not  concentrated  at  a single  point  or  even  in  a small  area  ol  the  cylinder 
but  extends  throughout  its  entire  width.  The  intensity  of  the  virtual  caustic  does  however  decrease  from  a 
maximum  on  the  \-axis  to  zero  on  the  y-axis.  In  some  cases,  c.g..  when  the  field  points  ol  interest  arc  tar 
from  the  cylinder  in  terms  of  cylinder  diameters,  it  will  be  possible  to  approximate  the  reflected  Held  by  a 
line  source  field  emanating  from  1 s(Figme  III  I). 

It  is  also  worth  noting  that  there  will  be  no  constructive  or  destructive  interference  between  the 
reflected  rays.  There  will,  of  course,  be  interfcience  between  the  direct  and  reflected  rays. 


IV.  ESTIMATE  OF  DIRECTLY  RADIATED  VERSUS  REFLECTED  POWER 


A long  circular  cylinder  irradiated  by  a plane  wave  at  right  angles  to  the  cylinder  axis  will  cause 
reflections  from  the  front  face  of  the  cylinder  according  to  the  well  known  equation  of  geometrical  optics 
(zero-order  approximation  of  C.TD): 


IC,I' 


>2 

a cos  0, 


-1/2 


(IV- 1 ) 


where  a ....  radius  of  the  cylinder 

r1  ...  distance  to  the  field  point 

0 angle  of  incidence 

\Ur\  reflected  field,  magnitude 

.10  . ..  amplitude  of  the  incident  field 

In  this  case,  both  constructive  and  destructive  interference  will  be  of  importance  depending  on  the  field 
point  ( F)  position.  If  we  call  /fi  the  absolute  value  of  the  ratio 

Idirect  field  I 
2 Irellected  field  I 


we  obtain  from  (IV-I  ) 


( IV- 2 > 


Note  that  for  r^/a  » I it  follows  that  /(r  » I . 

In  these  field  points  the  reflected  energy  may  be  neglected  compared  to  the  direct  radiation.  The 
value  of  IWjl  IS  shown  in  Table  IV-I  below: 


Table  IV  1 

(Direct  Fieldl 

lA  ■>  1 

2 IReflccted  Fieldl 

\rifa 

IX. 

1 

-> 

5 

10 

20 

50 

20 

1 769 

2.291 

1.412 

4.721 

6.600 

10.16 

40 

1 .900 

2.49  t 

1.749 

5.207 

7.295 

11.47 

60 

2.2.16 

1.000 

4.581 

6.401 

9.000 

14.18 

80 

1.518 

4.901 

7.654 

10.78 

15.21 

24.02 

Nolo  that  tor  large  values  of  </>  and/or  r-^  'a,  li 2 becomes  large  and  in  many  cases  the  reflected  field  may  he 
neglected  compared  to  the  direct  field.  Since  the  fields  will  add  vectorially  and  will  in  general  be  at  phase 
angles  different  from  0 or  180°.  the  effect  of  the  smaller  field  will  be  further  diminished.  The  extension  of 
this  analysis  to  oblique  incidence  and  noncircular,  but  convex  and  smooth  cylinders  is  possible  with  some 
further  work.  Also,  applying  the  plane  wave  spectrum  method,  nonplanar  incident  fields  mav  be  treated 
following  a similar  line  of  reasoning. 


V OBLIQUE  INCIDENCE.  REELECTED  POWER  VERSUS  DIRECT  RADI  \ HON 


In  the  interest  of  clarity,  the  proposed  approximations  have  been  presented  in  Sections  II  and  IV  for 
the  simplest  case  of  a plane  wave  at  perpendicular  incidence  on  a circular  cylinder.  We  have  found  that  in 
many  circumstances  when  moderate  accuracy  is  sufficient  very  simple  algebraic  expressions  w ill  permit  us  to 
decide  whether  the  diffraction  and  reflection  due  to  the  cylinder  may  be  neglected,  compared  to  other 
field  contributions,  in  a given  field  point.  The  above  cases  could  be  treated  with  two-dimensional  geometry. 

In  the  case  oi  a plane  wave  at  oblique  incidence  to  the  cylinder  axis,  we  shall  need  three-dimensional 
geometry  to  describe  the  problem.  Whereas  in  the  problem  of  Sections  II  and  IV  the  plane  of  incidence 
intersected  the  cylinder  in  a circle,  at  oblique  incidence  the  figure  of  intersection  will  be  an  ellipse.  The 
general  case  can  be  reduced  to  a special  (canonical)  case  if  the  coordinate  system  is  chosen  in  a special 
manner.  It  is  assumed  that  a contribution  of  the  field  at  the  field  point  F is  obtained  bv  reflection  from  the 
cylinder  t ailing  the  point  of  reflection  t).  the  \ -axis  of  the  chosen  coordinate  sy  stem  will  pass  through  Q 
while  the  /-axis  remains  unchanged,  fur  a given  direction  of  incidence  and  given  I the  point  0 can  be 
determined  An  appropriate  coordinate  transformation  will  then  produce  the  canonical  case.  The  reflected 
power  as  well  as  the  width  of  the  penumbra  region  can  then  be  determined  easily,  f rom  them  the  ratio  AS 
(Section  IV)  can  be  determined.  The  flowchart  for  this  procedure  is  given  in  Appendix  A2. 


VI  COMPARISON  WITH  EXPERIMENT 


I or  verification  of  the  physical  optics  method  used  in  Sections  IV  and  V.  measurements  by  C'.C.  H. 

TangI * * 4  were  used  Figure  VI  I shows  the  exact  solution  and  measured  points  for  the  backscattering  widths 

of  infinitely  long  perfectly  conducting  circular  cylinders  with  I -fields  parallel  to  the  axis  (normali/.ed  as 

)• 


Figure  VI- 1 . Backscattering  from  Circular  Cylinders  after  Tang2 * 4 

The  physical  optics  method  may  be  employed  in  order  to  calculate  the  values  of  the  backscattering  widths 
given  in  Figure  VI- 1 . 'Die  calculation  was  carried  through  in  Appendix  A7.  The  results  of  the  simple  physi- 
cal optics  method  agreed  in  this  example  with  the  exact  modal  solution  and  the  measured  values  to  within  3%. 
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APPENDIX  A 1 

ESTIMATE  OF  THE  TOTAL  DIFFRACTED  POWER  vs.  TOTAL  INCIDENT  POWER 


i 


*lwtei/  is  the  width  of  the  penumbra  region  on  the  mantle  of  the  cylinder  (I- ig.  A 1-1  )d  = [lu^lk)'1' 
I Kelerence  4 ) let  the  ratio  !)< H))  he  denoted  K 


R ~ R j 


la/ 


WJV 0 ^ 


(Alt) 


Thus,  we  obtain  Table  ( A I - 1 ) 

Table  All.  Ratio  of  Reflected  Power/Diffracted  Power 


ka 

K | = (te/2)'/3 

K.ldB) 

50 

2.924 

1 ’ ' 

4.66 

20 

2.154 

3.33 

10 

1.710 

2.33 

5 

1 .357 

1.33 

i 

1 .000 

0.00 

As  expected,  the  ratio  ol  reflected  power  to  diffracted  power  decreases  with  decreasing  diameter/X 
Reflected  power  equals  diffracted  power  when  ka  = 2.  The  total  power  /'  involved  in  the  scattering  process 

/ / A-  + /'/)  is  obtained  by  combining  Hquations  (Al  l ) and  (AI-2). 


P- 


(A  I -4) 


Normalized  to  unit  incident  power  and  substituting  for  (2/ak)'13  from  equation  (AI-3) 


( A 1 -5 ) 


Ihe  total  scattered  power  /'and  the  total  diffracted  power  are  contrasted  in  Table  AI-2  for  constant  k = I 
and  normalized  to  unit  incident  power 


Table  AI-2.  Scattered  Power  /’  and 
Diffracted  Power  Pp 


ka 

PIPa 

PdIPq 

1‘rII‘u 

I’nir R 

50 

134.2 

34.20 

100 

0.342 

20 

58.56 

18.56 

40 

0.4(i4 

10 

31.70 

1 1.70 

20 

0.585 

5 

17.37 

7.37 

10 

0.737 

8.- 

4.- 

4 

1 .000 

Note  that,  although  the  relative  importance  of  the  diffracted  versus  rellected  power  increases  as  ka  decreases 
its  absolute  magnitude  decreases  lor  constant  A and  given  incident  power  density. 
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APPENDIX  A2 

FLOWCHART  FOR  THE  COMPUTATION  OF  THE 
REFLECTED  FIELD  IN  THE  CASE  OF 
OBLIQUE  INCIDENCE  OF  A PLANE  WAVE  ON  A CIRCULAR  CYLINDER 


I 


t' 


c 

et  = edirect  + ereflected 


Figure  A 2- 1 . Illustration  to  Computer  Program  F lowchart  of  Figure  A 2-3. 

The  above  figure  shows  the  configuration  of  the  incident  ray  and  the  ray  obliquely  reflected  from  the 
cylinder.  The  /-axis  points  out  of  the  plane  of  the  paper.  The  position  of  the  vector  /■  = Ob'  is  given.  Note 
that  for  any  given  vector  the  coordinate  system  may  be  chosen  so  that  Q = (>(Vy . Ty . 0);  reflection  from  a 
point  Q (/(jty.Cy,  Zq  ) will  be  obtained  by  shifting  the  coordinate  system  upward  or  downward  along 
the  / axis  to  obtain  the  situation  shown  in  Figure  A2-I . Further  after  <t>  is  calculated  a further  coordinate 
x-lnrmatinn.  a rotation  about  0.  is  performed  to  bring  the  point  Q onto  the  x-axis.  The  canonical  position 
thus  produced  is  shown  in  Figure  A2-2.  For  all  field  points  where  reflected  rays  pass,  the  canonical  cases 
may  be  represented  as  in  Figure  A2-2.  with  only  0 varying.  The  computation  flowchart  is  shown  in  Figure 
A2-3. 


1 1 


Q’(-a,0,0) 


x 


Figure  A2-2.  Canonical  Case 
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Figure  A2-3.  Computation  of  the  Field 
Reflected  From  a Cylinder 
into  the  Field  Point  /■  for 
Obli<|tie  Incidence 
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APPENDIX  A3 

COMPARISON  OF  THE  PHYSICAL  OPTICS  RESULTS  WITH 
EXPERIMENT  AND  EXACT  MODAL  SOLUTIONS  IN 
BACKSCATTERING  FROM  A VERY  LONG  CIRCULAR  CYLINDER 


Figure  A3-2.  Illustration  to  the  Simplified  Analysis  of 
Backscattering  from  a Circular  Cylinder 


( A3- 1 ) 


o‘  is  the  echoing  width,  hy  definition 

2rr  (scattered  power/unit  tingle) 

incident  power/unit  length  (along  the  cylinder  diam  iter ) 

ill  = adO 
Jv 

— = cos  0 
Jl 

tlv  = d I cos  0 = a cos  OdO 


dP  dP 

incident  power/unit  length  — — = rrz 

dv  a cos  OdO 


(A3 -2) 


dP/dy  = constant  = c ...  power/unit  length  (along.v) 

The  power  c incident  along  unit  length  of y causes  power  /)'  scattered  front  unit  length  of  the  cylinder 
contour: 


, dp  dp  dy  _ 

n = — = — * — = C cos  0 
' Jl  dy  Jl 


The  power  dp  incident  on  the  infinitesimal  part  of  the  contour  dl  will  be  scattered  into  an  angle  2 dO 

c cos  0 dO 


Hence  the  scattered  power/unit  angle  = 


2 dO 


(A3-3I 


and  trom  (A3-I  MA3-2).  and  (A3-3) 


2ttc  cos  0 "a 


For  backscatter  0=0.  and  Kquation  (A34)  becomes 

och*air 


(A34) 


(A3- 5) 


To  compare  with  Figure  VI-1  we  form 


h 

k0a‘h/A:0a  = — = tr  (front  (A3-5)) 

From  Figure  A3- 1 we  read:  k^'Jk^a  — 3.07  for  the  experimental  values  as  well  as  the  exact  theoretical 
values.  Note  that  in  this  case  the  value  obtained  with  our  simple  physical  optics  model  gives  the  precise 
theoretical  and  measured  values  to  within  3%. 
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